Lattice QCD Calculation of Flavor Singlet Matrix 
Elements and TV- TV Scattering Lengths 



Yoshinobu Kuramashi 

National Laboratory for High Energy Physics (KEK) 
Tsukuba, Ibaraki 305, Japan 

Abstract 

We report on our recent study of flavor singlet matrix elements calculated 
with the wall source method without gauge fixing. Results are presented for the 
mass difference between 7/ and pseudoscalar octet mesons, the ir-N a term and 
the proton axial vector matrix elements. An exploratory calculation of the N- 
N scattering lengths is also discussed including a phenomenological study with 
one-boson exchange potentials. 

1 Introduction 

Numerical simulations of lattice QCD have been successfully applied to calculate a wide 
of variety of physically interesting observables. However, the class of observables treated 
to date are mostly restricted to the flavor non-singlet quantities and little is known for 
the flavor singlet ones. The reason is the technical difficulty of calculating amplitudes 
involving disconnected quark loops required in the latter case. We have found, however, 
that the method of wall source without gauge fixing, previously developed to deal with 
a similar computational difficulty with full hadron four-point functions]!]], allows to 
overcome the problem. Using this method we have recently calculated the rj' meson 
mass||, [|, the tc-N a term|4[] and the flavor singlet axial vector matrix elements of 
proton ||. The first aim of this article is to present a summary of our findings on these 
quantities. 

Calculation of hadron scattering lengths is a first step going beyond a study of static 
properties of hadrons based on lattice QCD. Of particular phenomenological interest 
is the nucleon-nucieon scattering lengths since, unlike the pion-pion and pion-nucleon 
cases, their values are not constrained by chiral symmetry. The second aim of this article 
is to present our study of the subject, including a one-boson exchange model calculation 
on the quark mass dependence of the scattering lengths |6|, [7[|. 

All of our calculations are carried out in quenched QCD at (3 = 5.7 using the Wilson 
quark action with the hopping parameters K = 0.150 — 0.168 on a 12 3 x 20, 16 3 x 20 or 
20 3 x 20 lattice. Gauge configurations are generated by the pseudo heat-bath method 
at 1000 sweep intervals. Errors in the physical quantities are estimated by a jackknife 
procedure. We correct quark field normalization by Jl — ?>K/AK C with K c = 0.1694 



and employ the tadpole- improved one- loop expressions for the renormalization factors 
of the quark bilinear operators with ajjg(l/a) = 0.2207 for the coupling constant. 

2 rf meson mass 

The large value of the flavor singlet rf meson mass in comparison to the masses of octet 
pseudoscalar mesons constitutes the well-known U(l) problem^. The current view on 
the resolution of this problem is based on the 1/N C argument |10||: the mass splitting 
m 2 , = ml, —ml is generated through iteration of virtual quark loops in the rf propagator, 
each of which gives a factor m\j {p 2 + m§), and the U(l) anomaly is supposed to give a 
large value to m Q . Our aim is to evaluate m directly in lattice QCD and examine the 
relation between m and the topological charge. 

In quenched QCD the rf propagator only consists of the one-quark- loop amplitude 
having a single pole l/(p 2 + mf) and the two-quark-loop amplitude with a double pole 
m o/(P 2 + m s) 2 - This means that the mass splitting mo can be extracted from the ratio 
of the two amplitudes defined by 



() _ toW(0))2-iocp r± { 



where each amplitude is projected onto the zero momentum state. 

In order to evaluate the two-loop amplitude, we calculate the quark propagator 
G(n,t) = J2( n ",t") G(n,t; n",t") obtained with unit source at every space-time site with- 
out gauge fixing. The expression J2 n Tr[G(n, 0)75] J2n' Tr [G* (n', 2)75] equals the two- 
loop amplitude aside from extra contributions of gauge-variant non-local terms. How- 
ever, these must cancel out in the ensemble average. This technique has the advantage 
that only a single quark matrix inversion is required for each gauge configuration, which 
enables us to achieve high statistics. 

Our calculations of the //-octet mass splitting are made with four values of the 
hopping parameter, K = 0.164, 0.165, 0.1665 on an L = 12 lattice and K = 0.168 on 
an L = 16 lattice. For each hopping parameter we use 300, 240, 240 and 43 gauge 
configurations, respectively, to evaluate one- and two-loop amplitudes. In the last case 
three exceptional configurations are excluded, for which we fail to find the solution of 
the quark matrix inversion. 

We extract the mass splitting mo fitting the data of R(t) to the linear form ([!]) 
using the value of m§ taken from a single exponential fitting of the pion propagator (see 
Fig. 2 of Ref. 0). The results for m multiplied by yN f = \fd> to obtain the physical 
value are plotted in Fig. as a function of the quark mass m q = (1/K — l/K c )/2 
using K c = 0.1694. We observe a linear increase of mo as m q decreases toward the 
chiral limit. Employing a linear extrapolation for m to obtain the value at m q = 0, 
we find m = 0.518(25) or m = 751(39)MeV in physical units with a -1 = 1.45(3)GeV 
determined from the p meson mass. 

The value of mo we obtained may be compared with the empirical mass difference 
mg = mi, -(Am 2 K + 3ml + m 2 ri )/8 w (866MeV) 2 , which is the naive estimation ignoring a 
small mixing between rj and rf . While a precise agreement is not obtained, both values 
are sufficiently close to conclude that the //-octet mass splitting can be understood 
within quenched QCD. 
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Figure 1: mo multiplied by yNf = v^3 as a function of m q = (1/K — l/K c )/2. Left 
ordinate and bottom abscissa are in lattice units, while right ordinate and top abscissa 
are in physical units using a -1 = 1.45GeV. Solid line is a linear fit which extrapolates 
to m = 0.518(25) ( 751(39)MeV ) at m q = 0. 



Let us next consider the topological charge dependence of the //-octet mass splitting 
mo- The topological charge on the lattice is defined as 

Q = 12 W^ ReTr (^ WfVHK ( 2 ) 

OZ7T n 

where U^{n) is the plaquette in the \w plane at site n. We calculate Q by the cooling 
procedure [PH employing 25 cooling sweeps to wash out ultraviolet fluctuations. 

To check the anomalous Ward identity relation ml = 2Njx/f 2 derived by Witten 
and Veneziano in the large N c limit 01, we calculate the topological susceptibility x — 
(Q 2 ) I (space-time volume) using our 300 pure gauge configurations on an L = 12 lattice, 
and find x = 4.76(39) x 10~ 4 . Combining with f n = 0.0676(24) at m q = obtained on 
the same set of gauge configurations used for the calculation of m , we find yj 2Nfx/ f% = 
1146(67)MeV for N f = 3 with the aid of a' 1 = 1.45(3)GeV. The result is somewhat 
higher compared to the value of estimated from the direct measurement. It is possible 
that the disagreement arises from chiral symmetry breaking of the Wilson quark action 
which induces extra terms in the anomalous Ward identity. 

In order to further examine the relation between the //-octet mass splitting and the 
topological charge we classify gauge configurations according to the value of which 
is the absolute magnitude of Q rounded off to the nearest integer, and evaluate R(t; \Q\) 
on each ensemble. For the one-loop amplitude in the denominator we use the average 
over the whole ensemble since it shows little variation depending on \Q\. The mass 
splitting rriQ is extracted by fitting the data of R(t; \Q\) to the form of ([!]). 
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Figure 2: Representative values of the ratio R(t; \Q\) of two- and single-quark loop 
contributions to the rj propagator for several values of \Q\ at K — 0.1665 on an L = 12 
lattice. 



In FigJ| we present R(t; \Q\) at K = 0.1665 for typical values of \Q\: \Q\ = 0, 2, 5. 
We see that the function R(t; \Q\) is consistent with zero for \Q\ — and clearly increases 
as \Q\ increases. This fact demonstrates that the excitation of instantons generates the 
contribution of the two-loop amplitude. The dependence of the //-octet mass splitting 
m on \Q\ is listed in Table |1|. The increase of m with \Q\ shows that the rj meson 
mass increases with the winding numbers of the instantons. 

In full QCD with dynamical quarks topological excitations become suppressed toward 
vanishing sea quark mass due to the fermionic zero modes. It is an interesting problem 
to understand how the rj meson mass survives the suppression. Full QCD simulations 
for examining this point have to overcome the difficulty that topological fluctuations are 
known to exhibit long-range correlations in Monte Carlo time|]l 2||. 



Table 1: \Q\ dependence of tuq in lattice units at K = 0.1665 on an L = 12 lattice. 
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3 Flavor singlet nucleon matrix elements 
3.1 scalar matrix elements 

The 7C-N a term is a fundamental quantity measuring the degree of chiral symmetry 
breaking in the strong interaction and defined as a = (m u + m^)/2 • (N\uu + dd\N). 
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Table 2: Nucleoli scalar matrix elements corrected by tadpole-improved Z factor as a 
function of K for up and down quarks. (N\ss\N) represents an interpolated value at 
the physical strange quark mass K s = 0.1648 with up and down quark masses fixed. 



K 
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\ iV l**l iV / 
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conn. 


disc. 








0.1600 


12 


300 


2.323(15) 


3.56(76) 


2.09(25) 


0.749(5) 


-0.074(2) 




16 


260 


2.326(17) 


3.02(98) 


1.78(32) 


0.749(7) 


-0.075(6) 


0.1640 


12 


300 


2.413(30) 


4.58(92) 


2.36(29) 


0.762(10) 


-0.126(6) 




16 


260 


2.378(30) 


4.1(1.2) 


2.10(40) 


0.748(12) 


-0.128(11) 


0.1665 


12 


400 


2.693(81) 


5.1(1.1) 


2.66(35) 


0.831(25) 


-0.200(16) 




16 


260 


2.565(73) 


5.6(1.6) 


2.70(51) 


0.783(26) 


-0.204(21) 


K c 


12 




2.615(61) 


5.8(1.4) 


2.84(44) 


0.802(19) 


-0.208(11) 




16 




2.515(60) 


6.2(1.9) 


2.89(61) 


0.765(23) 


-0.222(20) 



Study of the a term has a long history and still shows some controversy due to the 
difficulty of estimating the sea quark contribution including the strangeness content of 
the nucleon. A direct calculation of the a term in lattice QCD potentially provides a 
thorough resolution of this issue. Some previous calculations O, |TJj suggest that the 



disconnected contribution is as large as the connected one. However, those calculations 
were made only at quite heavy quark masses and it seems difficult to extrapolate their 
result to the physical point (see Fig. 1 of Ref. []i~5]| ). 

We have made a systematic calculation of the quark mass dependence of the 7r- 
N a term. The nucleon matrix elements of scalar quark density are extracted from 
a linear fit in t of the ratio of nucleon three-point function divided by the nucleon 
propagator [16]. Dirichlet boundary condition is employed in the temporal direction for 
quark propagators. The source method [|T7f is used for the connected contribution. The 
disconnected quark loops are evaluated using quark propagators solved with unit source 
for every space time site without gauge fixing. 

In Table ^| we summarize our results for the nucleon scalar matrix elements for each 
spatial size L and the hopping parameter K. Since the results for L = 12 and 16 
shows little finite size effects, we give the averages of L = 12 and L = 16 results in the 
following. 

For the isoscalar matrix element (N\uu + dd\N) we observe that the disconnected 
contribution is about twice larger than the connected one. In the the chiral limit we 
find for the ratio cr disc /cr conn = 2.35(46). 

Calculating the physical value of the a term requires an estimate of the quark mass. 
For this purpose hadron masses m n , m p , tun are calculated using the same set of hopping 
parameters on the same set of gauge configurations. The light quark mass m = (m u + 
md)/2 = 0.0034(1) (m q = (1/K — l/K c )/2) is determined using m n /m p = 0.18 and the 
strange quark mass m s = 0.0829(19) (K s = 0.1648) from m^jm p = 0.64. The lattice 
scale a -1 = 1.46(2)GeV is fixed by m p at K c . Here we take the average of results for the 
L — 12 and 16 lattices since we do not detect any finite size effects for hadron masses 
between the two lattice sizes fi~5|j. 

We estimate the physical value of a in two ways. The isoscalar matrix element at 
K c multiplied by the light quark mass rh gives a = 44(6)MeV. Alternatively, the ratio 
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m^a/m^ extrapolated linearly in m q to the chiral limit yields 60(9)MeV with the aid 
of experimental pion and nucleon masses. The difference between the two values of 
a originates from the large ratio vn^jmp = 1.47(3) in our simulation carried out at 
/3 — 5.7. Within this uncertainty these estimates show an encouraging agreement with 
experimental estimates a = 45MeV]TB[. For the strangeness content of nucleon we find 
a fairly large value for the y parameter y = 2(N\ss\N) / (N\uu + dd\N) = 0.66(13), as 
compared to a phenomeno logical estimate y ~ 0.2 JT/5| . 

Obvious systematic uncertainties in our results are scaling violation and sea quark 
effects. In addition we should note the potential problem that the physical light quark 
mass for the Wilson quark action substantially decreases toward smaller lattice spacing 
in quenched QCD and that the full QCD values are factor 2 — 3 smaller than those for 
quenched QCDllJ. 

Let us also make a remark concerning the flavor non-singlet SU (3) reduced matrix el- 
ements Fs and D$- These may be evaluated from the relations Fs = (p\uu — ss\p) /2 and 
Ds = (p\uu — 2dd + ss\p) /2 with \p) the proton state, where we may ignore disconnected 
contributions to leading order in flavor symmetry breaking. For the mass splitting of 
the baryon octet, our results for the reduced matrix elements yield m~ — ttin = {ni s — 
rh) ■ 2F S = 0.182(6)GeV and -3/2(m E - m A ) = (m s - m) • 2D s =-0.0498(30)GeV. The 
magnitude of these values is roughly a factor two smaller compared to the experimental 
values 0.379GeV and -0.116GeV, while our result for the ratio D s /F s = -0.275(16) 
is reasonably consistent with the experimental estimate Ds/Fs ~ —0.32. A source of 
uncertainty, besides scaling violations and sea quark effects, is the validity of the first 
order approximation for the mass splitting. This point should be resolved by a simul- 
taneous calculation of the baryon mass splitting and scalar matrix elements, which we 
leave for future analysis. 



3.2 axial vector matrix elements 

Interest in the flavor singlet axial vector matrix elements of proton was inspired by the 
EMC result for the spin dependent structure function gi of proton [^0J, which indicated 
that the fraction of proton spin carried by quarks has a small value and that the strange 
quark contribution is unexpectedly large and negative. Recent analyses combining data 



Table 3: Proton axial vector matrix elements corrected by tadpole- improved Z factor 
as a function of K for up and down quarks. Ad disc equals Aw disc . As represents an 
interpolated value at the physical strange quark mass K s = 0.1648 with up and down 
quark masses fixed. 
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0.3974(55) 0.6358(83) 




16 
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-0.066(23) 
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0.1665 


12 


0.761(20) 


-0.232(11) 






0.3806(98) 0.611(15) 




16 


0.818(39) 


-0.231(23) 


-0.093(54) 


-0.087(33) 


0.409(20) 0.641(16) 
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12 


0.694(20) 


-0.2284(96) 






0.347(10) 0.575(15) 




16 
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Figure 3: Axial coupling g& for quenched (open) and full QCD (closed) with Wilson 
quark action. Our results for L = 16 (circles) and L = 12 (up-triangles) are at (3 = 5.7. 
(3 = 6.0 for other quenched results (open squares |T3 
Full QCD results (closed squares |T3| 



diamonds 29 



are at (3 

the upper and lower bound of la for the experimental value gA 



down-triangles [| 
5.4 — 5.6. Horizontal lines represent 
1.2573 ±0.0028. 



IT], [Z2JJ , deuteron[T2^, [23[] and neutron |f24[| targets yield AS = Au + 
-0.83(3) - 0.43(3) - 0.10(3) = 0.31(7) at the renormalization point /i 2 = 

0.31(4) at fi 2 = oof26j, where Aq is 



on proton pO 
Ad + As = 

10GeV 2 [[25| or +0.832(15) - 0.425(15) - 0.097(18) = u.6L{4) at yr = ooj 
the forward matrix element of cpiCi^q for a proton polarized in the i-th direction averaged 
over % = 1, 2, 3. These results are far from the naive constituent quark model which 
predict Au = 4/3, Ad = -1/3 and As = 0. 

We have attempted a lattice calculation of Ag for each flavor q including both con- 
nected and disconnected contributions. We evaluate the proton matrix elements of the 
axial vector current employing the same calculational technique, the same set of gauge 
configurations and the same hopping parameters as for the case of the nucleon scalar 
matrix elements. 

Our results are summarized in Table |3|. We observe little finite size effects for the 
connected contributions. For the disconnected piece statistical fluctuations turned out 
to be significantly larger than those for the scalar matrix elements, to the extent that a 
reliable signal was not obtained on an L = 12 lattice. 

Our results for L = 16 show that the disconnected contribution, albeit involving 
a substantial error of 50%, are negative and that their magnitude is small compared 
to those for the connected ones. Also, the disconnected contributions slightly increase 
toward the chiral limit while the connected ones exhibit a slight decrease. Adding both 
contributions for each flavor and making a linear extrapolation to K — K c , we find 
AS = An + Ad + As = +0.638(54) - 0.347(46) - 0.109(30) = +0.18(10) with the 
disconnected contribution to Am and Ad equal to —0.119(44). These values show a 



reasonable agreement with the phenomenological estimates quoted above[25, 26 
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In Table |3| we also list the flavor non-singlet matrix elements Fa = (Aw — As)/2 and 
Da = (Ait — 2Ag?+As)/2 evaluated without the disconnected contributions under an as- 
sumption of exact SU(3) flavor symmetry. In the chiral limit we find Fa = 0.382(18) and 
D A = 0.607(14) for L = 16, which give the axial charge Qa — Fa + D A = 0.985(25) and 
the ratio Fa/Da = 0.629(33). Compared to the experimental values g A = 1.2573(28) 
and Fa/ Da = 0.586(19) |27|], the ratio shows a good agreement, while g A is about 25% 
smaller. Other calculations of for quenched QCD at (5 = 6.0 |TB], |23|, |2P| and for full 



QCD[|13| at (3 = 5.4 — 5.6, with the lattice spacing of a ~ 0.15 — O.lfm, yield similar 
results if we use the same renormalization factor as we employed. We illustrate this 
point in Fig. |3|. The agreement of results for a variety of simulation parameters indi- 
cates that possible systematic errors such as finite size effects, scaling violations and 
quenching effects are too small to explain the deviation of the lattice estimates from the 
experimental value. The origin of the discrepancy is an open problem which should be 
resolved in future work. 



4 N-N scattering lengths 

Lattice calculation of hadron scattering lengths is based on the Liischer's formula pOf 
which relates the s-wave scattering lengths ao to the energy shift of two particle state 
at zero relative momentum confined in a finite periodic spatial box of a size L 3 : 

77i / , \ 27r(mi + m 2 )a / / a \ f a o\ 2 \ . nlT -^ 

with c\ = —2.837297, c% = 6.375183. For the nucleon-nucleon case application of this 
formula faces the difficulty that the experimental values ao( 3 5i) = — 5.432(5)fm and 
a o( 1 5o) — +20.1(4)fm|l3T| are quite large so that a large lattice is required to suppress 



0(L~ 6 ) corrections in (^). The negative sign in the 3 Si channel due to the presence of 
the deuteron bound state (Levinson's theorem) brings an additional complication that 
the energy of the lowest scattering state orthogonal to the bound-state deuteron has to 
be computed to apply (|3p. 

To avoid this complication we take a strategy to start our simulation from the heavy 
quark mass region. In order to have an idea on the behavior of N-N interactions toward 
large quark masses, we carry out a phenomenological analysis employing a model of one- 
boson exchange potentials. Since these models include a large number of parameters 
fine-tuned to reproduce the experimental results for the scattering phase shifts, we make 
this study not to analyze detailed quark mass dependence of scattering lengths, but to 
obtain a rough feature of the nuclear force under variation of quark mass around the 
physical one. 

Specifically we employ the model of Ref. |Q and vary the value of according to 



our lattice results for the slope m^/m q . The mass of the scalar particle a is also varied 
in proportion to m n since a represents the contribution of two-pion exchange. For N, 
p and uj masses we assume the linear m q dependence m h = a>h + bh ■ m q , fixing the 
slope bh with the aid of our lattice results. Other meson masses m v , ms, and the 
meson-nucleon coupling constants are fixed to the the physical value. 

The quark mass dependence of the N-N scattering lengths obtained in this way is 
shown in Fig. HI. We find a divergence of the scattering length in the 3 Si channel taking 
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place at m q = 6.3MeV that signals unbinding of the deuteron, which is only 30% larger 
than the physical point m q = 4.8MeV. On the other hand the 1 5'o channel forms the 
bound state at a slightly smaller quark mass of m q = 4.6MeV. For both channels the 
magnitude of the scattering length reduces quickly for heavier quark and becomes close 
to the order of 1 fm. It is reasonable to expect that the scattering lengths take values 
similar to the hadron size in the absence of bound-state effects. 

These results suggest that the deuteron becomes unbound as the quark mass in- 
creases, and hence Luscher's formula can be applied to extract the N-N scattering 
lengths for heavy quark. Our lattice study of the N-N scattering lengths is carried out 
based on this consideration using the same technique as we employed for the tt-tt and tt- 
N cases P, 0]. For the calculation of the nucleon four-point function gauge configurations 
are fixed to the Coulomb gauge over all space-time to enhance signals. We analyzed 
20, 30 and 20 gauge configurations for the hopping parameters K = 0.150, 0.160 and 
0.164, respectively, on a 20 3 x 20 lattice. In Fig. [5] we compare the results for the N-N 
scattering lengths with those for tt-tt and tt-N cases. As expected, the N-N scattering 
lengths are of the order of 1 fm with a positive sign for both eigenchannels. They are 
substantially larger than the tt-tt and tt-N scattering lengths already for a heavy quark 
corresponding to m^/nip m 0.74. Also we note the trend, albeit with sizable errors, that 
the values for the 3 Si channel are larger than those for the 1 So channel. This indicates 
a stronger attraction in the 3 Si channel, consistent with the formation of the deuteron 
bound state at the physical quark mass where a bound state is not yet formed in the 
1 So channel. 

If we reduce the quark mass toward the physical point, we expect an increase of the 
N-N scattering lengths in both eigenchannels. For the 3 Si channel the scattering length 
should diverge at some quark mass signalling that the deuteron bound state is formed. 
However, actually observing such a phenomena will be difficult. Large lattice sizes far 
exceeding the size of 2 — 3fm, which are accessible in current lattice QCD simulations, 
will be needed. Furthermore increase of statistical fluctuations for large times in the 
N-N four-point functions is expected to become quite severe toward the chiral limit as 
shown in Fig. 23 of Ref. M. 

5 Conclusions 

The method of wall source without gauge fixing has proven to be an effective tool for 
exploring the role of disconnected quark loop amplitudes in the flavor singlet sector of 
strong interactions. Our results for the r/'-octet mass splitting demonstrate that topo- 
logically non-trivial gauge configurations boost the r/ meson mass to the level observed 
experimentally. Our application of the method to flavor singlet nucleon matrix elements 
have shown that disconnected contributions are substantial for the tt-N sigma term and 
the quark content of proton spin. The total contribution we found for these quanti- 
ties are reasonably consistent with the experimental estimates. On the other hand, the 
flavor non-singlet pieces calculated from the connected contributions show substantial 
discrepancies between the lattice estimates and experiments. 

We have also described our exploratory calculation of the N-N scattering lengths at 
heavy quark masses, for which we found results expected from considerations based on 
a phenomenological model of one-boson exchange potentials. 
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